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Abstract 



o 

<D 

^0 " We study diffusion and wave equations in networks. Combining semigroup and 

variational methods we obtain well-posedness and many nice properties of the solutions 
in general L p -context. Following earlier articles of other authors, we discuss how the 
spectrum of the generator can be connected to the structure of the network. We 
conclude by describing asymptotic behavior of solutions to the diffusion problem. 

^ ! 1 Introduction 

In this paper we continue the study of dynamical processes in networks using semigroup 
methods. While [KS05| and |MS06j studied flow and transport processes, the aim of the 
present paper is to combine variational and semigroup methods in order to obtain the well- 
posedness of initial value problems associated with diffusion and wave equations. We thus 
Ijr^ ' consider first and second order problems 

OV 

iij(t,x) = {cju'j)' (£, x) and ilj(t,x) = [cju'^)' [t,x), t > 0, x € (0,1), 



. where cy(-) and Uj(t,-) are functions on parameterized edges e^ of a finite network. The 

node conditions in (|2.ip and (|3.1[) below impose continuity and Kirchhoff laws in the ram- 
C\| . ification vertices. Problems of this kind have already been treated by many authors both 

from the mathematical and physical communities - among others, we mention the earlier 
articles of Lumer }Lu80| . Ali Mehmeti [ATM] . Roth |Ro84j . von Below |Be85j . Nicaise [Ni85] . 
Exner |Ex89| . Cattaneo [Ca97], Kottos and Smilansky |KS97j . Kostrykin and Schrader jKS99j . 
and Kuchment jKu02j . as well as the monographs |Ni93] . |A194j . and |LLS94| . and the pro- 
ceedings jABNOlj . 

Since the pioneering work of Beurling and Deny in the 1950s, variational methods have 
been greatly developed. In combination with the theory of strongly continuous semigroups of 
operators, they provide a powerful tool to discuss properties of solutions to many parabolic 
and hyperbolic problems, cf. [Da89] . |Ar04] . and |Ou04j . While L 2 -techniques like the lemma 
of Lax-Milgram have been used in most of the above mentioned papers, our paper seems to 
be the first applying variational methods to obtain positivity, ultracontractivity, and stability 
for network equations in a general L p -context. This is the main aim of Sections 2 and 3. We 
remark that positivity of the semigroup governing the diffusion problem with much more 
general nodal conditions has been characterized, by algebraic methods, in |KS06j . 

We then proceed to study the qualitative behavior of the solutions. To that purpose 
we obtain in Section 4 a characteristic equation for the spectrum of the generator and de- 
scribe the appropriate eigensolutions. We reprove some results from |Be85| . |Ni85| . |Ni87bj . 
and |Be88bj in our setting with slight generalizations. We see that the spectrum is deter- 
mined by the structure of the network and corresponds to the spectrum of the Laplacian 
matrix known from graph theory, see |Mo91j . We give an explicit connection between the 
two spectra and show the impact of this to our problem. This relates to the well-known 
question: "Can one hear the shape of a drum?" , first addressed by Kac in |Ka66l . Concern- 
ing differential operators on graphs, the analogous question "Can one her the shape of a 
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network?" has been formulated and answered in the negative by von Below in his contribu- 
tion to [ABNOlj . Quite surprisingly, the same question was raised at the same time in the 
almost homonymous paper [GSOlj by Guttkin and Smilansky. They answered it in the pos- 
itive, by studying Schrodingcr operator on a finite, simple graph with rationally independent 
arc lengths and imposing some further assumptions on matching conditions at the vertices. 
In graph theory, however, it is well known that there are many graphs sharing the same 
spectrum, see |DH03j . Also in our case, the spectrum itself does not determine the network. 

In Section 5 we study the asymptotic behavior of solutions to the diffusion problem. To 
our knowledge this topic has not yet been properly treated by other authors. We show that 
the solutions always converge towards an equilibrium with rate of convergence depending 
on the structure of the network. This is discussed for special classes of networks. In similar 
contexts, convergence to equilibria has already been discussed, e.g., in |BN96j . 



2 The wave equation on a network 

We consider a finite connected network, represented by a finite graph G with m edges 
ei, . . . , e m and n vertices vi, . . . , v„. We assume that all the vertices have degree at least 
2, i.e., that each vertex is incident to at least 2 edges. Furthermore, we assume that G is 
simple, that is it has no multiple edges or loops. We normalize and parameterize the edges 
on the interval [0,1]. The structure of the network is given by the nxm matrices $ + :— (cfyfj ) 
and $ _ := (<frij) defined by 

b + . J 1, if e, (! 



0, otherwise, anC ^ Tl i ' 1 0, otherwise. 



We refer to [KS0 5 for terminology. The nxm matrix $ := (4>ij) defined by 

$:=$+- $- 

is known in graph theory as incidence matrix of the graph G. Further, let r(v^) be the set 
of all the indices of the edges having an endpoint at v,;, i.e., 

r(vi) := {j G {1, . . .,m} : e^O) = or e,-(l) = vj . 

For the sake of simplicity, we denote the value of the functions Cj(-) and Uj{t, •) at or 1 by 
Cj(vi) and Uj(t,Vi), if ej(0) = v.; or ej(l) = v^, respectively. With an abuse of notation, we 
also set u'j(t,Vi) = Cj(\/i) := whenever j ^ T(vi). When convenient, we shall also write the 
functions Uj in vector form, i.e., u = (ui, . . . ,u m ) T . 

Wc start with the second order problem 
(2.1) 

Uj{t,x) = (cju'j) 1 '(t,x), t G M, x G (0, 1), j = l,...,m, (a) 
Uj(t,Vi) = u e (t,Vi), t G M, j,£ G r(vi), i = l,...,n, (b) 
YJjLi&oViCjiviju'^t^i) = 0, t G M, i = l,...,n, (c) 

tij(0, a;) = fj{x), x G (0, 1), j = 1, . . . , m, (d) 

Uj(0,x) = gj(x), x G (0,1), j = l,...,m, (e) 

on the network. Note that Cj(-) and Uj(t, ■) are functions on the edge e^ of the network, so 
that the right-hand side of (12.1b ) reads in fact as 

(cju'jj'it, •) = A (c^u^j (t, ■), t G R, j = 1, . . . ,m. 

The functions c\, . . . , c m are the weights of the edges, and throughout this section wc 
assume that < Cj G i/ 1 (0,l), j = l,...,m. They represent the different speeds of 
propagation along each edge of the network G. The equation (12.1b ) represents the continuity 
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of the values attained by the system at the vertices. The coefficients jij, j = 1, ...,m, 
are strictly positive constants that influence the distribution of impulse happening in the 
ramification nodes according to the Kirchhoff-typc law (|2.f b ). 

We now introduce weighted incidence matrices := (wj) and <I>~ := (<^-) with entries 

w + . = { ft c j>'.)> ife 3 (0)=vi, w - . = / WCjK), ifej(l) = Vi, 

y ' ' [ 0, otherwise, 13 ' 1 0, otherwise. 

With these notations, equation (|2.f b ) can be rewritten as 

(2.2) 3d g C" s.t. ($ + ) T d = u(i, 0) and ($~) T d = u(i, 1 ), f g K, 

while the Kirchhoff law (|2.1b ) becomes 

*+«'(*, o) = *-«'(*, i)> ten. 

We are now in the position to rewrite our system in form of a second order abstract 
Cauchy problem. First we consider the (complex) Hilbert space 

m 

X 2 :=Y\L 2 {Q,l-^ J dx) 



endowed with the natural inner product 



(f,9)x 2 -=Y1 / fj( x )9j(x)fijdx, 
.7=1 J ° 



Observe thatX2 is isomorphic to (L 2 (0, l)) m with equivalence of norms. Moreover, X 2 
is in fact a Hilbert lattice whose positive cone consists of m copies of the positive cone of 
L 2 (0, 1; /jjdx) sa L 2 (0, 1). On X 2 we define an operator 

(2.3) A:= 

with domain 

(2 4) := {/ g (i? 2 (0, 1))™ : $+/'(0) = and 

3d g C" s.t. ($+) T rf = /(0) and ($-) T d = /(f)} . 

With this notations, we can finally rewrite (|2.f j) in form of a second order abstract 
Cauchy problem 

= Au(i), t g K, 
«(0) = f, 
6(0) = g, 

on X 2 . By means of variational techniques, we are going to show that A enjoys several nice 
properties. We follow the techniques of |Da89j and |ABHN011 Sec. 7.1]. 

Lemma 2.1. Consider the sesquilinear form 



3=1 

on the Hilbert space X 2 with domain 
D 



m -l 



(a) =V:=(f€ (H\0, l)) m : 3d g C" s.t. ($+) T d = /(0) and ($") T d = /(l)} 
TTien a is densely defined and has the following properties: 
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• (symmetry) : a(f,g) = a(g, f) for all f,g£ D(a), 

• (positivity) : a(f, /) > for all f g D(a), 

• (closcdncss) : V is complete for the form norm \\f\\ a := <\J a (fi f) + ll/ll x 2 > 

• (continuity) : \a(f,g)\ < -W||/||a||fl , ||a for some M > and all f,g g D(a). 

Proof. It is apparent that V is a linear subspacc of X2- Observe that (C£°(0, 1))™ C V . 
It follows that V is dense in X2, as by definition L 2 (0, 1) is the closure of C£°(0, 1) in the 
L 2 -norm. By assumption, the weights Cj are strictly positive, so that in particular o is 
symmetric and also positive, since 



o(/,/) = E / wiWiix^dx 

.7 = 1 ^ 



> for all / g V. 



Furthermore, 1/ becomes a Hilbert space whenever equipped with the inner product 

m -1 

(/.S)v := X] / + fj{ x WP)) Hjdx, f,g g V, 

since V is a closed subspace of (ff 1 (0, 1))™. Set 

c := min min cAx), C := max max c,(x). 

l<j<mie[0,l] l<i<mie[0,l] 

Then one has 

(cA^ll/ll^ll/H^^vi)!!/!! 2 ,, fev, 

so that the form norm || • ||„ is equivalent to the norm j| • \\y. Since V is complete with 
respect to || • \\v, the closedness of a follows at once. 

Finally, a is continuous. To see this, take f,g £ V and observe that 



m „1 

|a(/, fl )| < CJ2 \ / ^r j {x)g' j {x)dx\ 
3=1 Jo 

m 

< ll/jlU 2 (0,l;Mi*B)llfl r j'IU 2 (0,i;w<te) 



c I m \ 2 I m 

— "2 I ll/jlli 3 ro,l;m<fa^ J j Hflilli^O.lijUj-dg) 
2 • (c A 1) 

by the Cauchy-Schwartz inequality. □ 

Definition 2.2. From the form a we can obtain a unique operator (B, D(B)) in the following 
way: 

D{B) := {/ g V : 3g g X 2 s.i. a(/, ft) = ( 5 , ft)x 2 V/i g F} , 
5/ := -.9- 

FFe say that the operator (B, D(B)) is associated with the form a. 

Lemma 2.3. The operator associated with the form a is (A,D(A)) defined in (|2.3j) - (|2.4|) . 
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Proof. Denote by (B, D{B)) the operator associated with a. Let us first show that Ad B. 
Take / G D(A). Then for all ft G V 



*(/> h) = Ejli Jo 1 WCj (a;)^- (x)dx 



(2.6) 



= Ejli [MjCj/j7v] - £™ i Jo 1 M( c jfj)'( x ) h j( x ) dx - 
Using the incidence matrix $ = $ + — $ _ , the first term above can be written as 

rn m n 

X! [/ i J c i-/j^] I!X''''i iv ^ ., _ <t>ij)fj(vi)hj{\>i). 



3=1 



j = l i=l 



Observe now that the condition 

3d G C n s.t. ($+) T d = ft(0) and ($~) T <i = ft(l) 

in the definition of V implies that ft is continuous in the vertices, i.e., there exist d±, . . . ,d n G 
C such that ftj(vj) = di for all j G r(vj), i = 1, . . . ,n. Summing up and using the other 
condition $+/'(0) = in 13(A) we obtain that 

m 

= E?=i 5X>« - <4)/j(vi) - E™i lii^yixjh^jujdx 

3=1 

= -(A/, ft)x 2 , 

which makes sense because A/ G X 2 . The proof of the inclusion A C B is completed. 

To check the converse inclusion B d A take / G D(B). By definition, there exists g G X2 
such that 

(2-7) X/ (J, j c j (x)f' j (x)h' j (x)dx = a(f,h) = {g,h)x 2 =y2 9j(x)hj(x)fj,jdx 
„-_i Jo „•_-, Jo 



3=1 



3=1 ' 



for all ft G V, hence in particular for all ft-? G V of the form 

j th row, G ffo(0,l). 



h 

From this follows that (|2.7II in fact implies 



fXjCj(x)f'j(x)h'Ax)dx= / gj(x)hj(x)fijdx for all j = 1, . . . , m, hj £ Hq(0,1). 
Jo 

By definition of weak derivative this means that Cj ■ /j G H 1 (Q, 1) for all j = 1, . . . , m. Since 
< Cj G -ff 1 ^ 1)) ^ follows that in fact /j G -ff 1 ^ 1) for all j = 1, . . . , m. We conclude 
that / G (i? 2 (0, 1)) . Moreover, integrating by parts as in (|2.6|) we see that if (|2.7|) holds 
for some ft G V, then 



71 771 



1=1 j=i 

where is the joint value attained at the vertex by all hj, j G r(vj). Since ft 6 7 is 
arbitrary, this means that 

XKj ~ w «)/j( v = for alH = 1, . . . , n, 

3=1 
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that is, $+/'(0) = $~/'(l). Therefore / € -D(-A) and 

"X!/ N( c ifi)'{ x ) h >i( x ) dx = Y) / 9j( x )hj(x)fj,jdx 

3 = 1 J ° 3 = 1 J ° 

holds for all /i £ V. This implies that A/ = —g, and the proof is complete. □ 

We are now able to use some well-known results on sesquilinear forms (cf. [ABHNOll 
iDaMIUuOlj ) in order to obtain nice properties of our operator A. 

Proposition 2.4. The operator (A,D(A)) defined in (|2 . 3[) — (|2 .4p is self-adjoint and dissi- 
pative. Thus, it generates a cosine operator function with associated phase space V x X2. 

Proof. By Lemmas 12.11 and 12.31 we are in the situation described in [ABHN01, Sec. 7.1] 
for H = X2, V = D(a), (-\-)v = •); <*> = 1, and A = Ah- Thus the claim follows 
by [ABHNOll Proposition 7.1.1], [ABHNOll Example 7.1.2] and the fact that self-adjoint 
operators are unitarily equivalent to multiplication operators (see also the remark at p. 413 
in |ABHJN01] n ). □ 

We can now state the main result of this section. This generalizes the well-poscdncss 
and regularity results in |A184] , [A186J , [Be88[ , and |CF03j , where only the case of constant 
or smooth coefficients c\ , . . . , c m was considered 



Theorem 2.5. The problem (|2.1I) is well-posed, i.e., for all f € V and g G X2 it admits a 
unique mild solution that continuously depends on the initial data. 

If further cj G C°°[0, 1], j = 1, . . . ,m, and the initial conditions f,g £ (C£°[0, l]) m , then 
the solution is of class (C°°[0, 1])™. 

Proof. It is well-known (see e.g. [ABHNOll Cor. 3.14.12]) that 

u(t):=C{t,A)f + S(t,A)g, teR, 

yields the unique mild solution to ([2.5[1 for all initial data (f , g) in the phase space, where 
we denote by (C(t, A))teu and (S(t, A)) t ^m the cosine and sine operator functions generated 
by A, respectively. The assertion about regularity of solutions follows directly from basic 
properties of cosine and sine operator functions. □ 



3 The heat equation on a network 

We now consider again the same network G and, under the same assumptions and with the 
same notations of Section 2, we turn our attention to the first order problem 

(cjuj.)' *>0, xG (0,1), j = l,...,m, 

u e (t,Vi), t>0, j,i€T(vi), i = l,...,n, 

0, t > 0, i = 1, . . . , n, 

fj(x), x G (0,1), j = l,...,m, 

This equation describes a diffusion process that takes place in a network and c\ , . . . , c m G 
C 1 [0, 1] are (variable) diffusion coefficients or conductances. Again, we are imposing conti- 
nuity and Kirchhoff-type conditions in the ramification nodes (controlled by some constants 

Mil • • • j Mm)- 

It is already known that such a problem is well-posed in an i 2 -context, cf. |Be88j . 
Moreover, at least for the case of constant weights c\, . . . , c m and [i\ = . . . = \i m = 1 the 
heat kernel has been computed in (Ni87^, thus yielding well-posedness in other L p -spaces. 
We show by variational methods that the semigroup governing (j3.1[) is L°°-contractive, 
and hence we can extend the well-posedness result to an L p -context by interpolation in 



(3.1) 



iij(t, x) 

Uj{t,Mi) 

Uj(0,x) 



G 



the general case of variable diffusion coefficients. In particular, the analyticity of the L p - 
semigroups seems to be a new result. Also observe that, by the bounded perturbation 
theorem, this also yields well-poscdncss for the Cauchy problem associated to the analogous 
cable equation, cf. |Ni87b] . 
Let 

m 

X p := \\_LP{Q,l;njdx), p€[l,oo]. 

3=1 

We have already seen in Proposition 12.41 that A is a self-adjoint and dissipative operator on 
X%. By the spectral theorem, this shows that A generates a contractive, analytic semigroup 
of angle ^ , and in particular the first order abstract Cauchy problem 

J u{t) = Au{t), t > 0, 
1 u(0) = f, 

is well-posed in X 2 . In fact, much more can be said. 

Lemma 3.1. The semigroup (l2(^))t>o on X 2 , associated with a, is sub-Markovian, i.e., it 
is real, positive, and contractive on X^ . 

Proof. By |Ou04[ Prop. 2.5, Thm. 2.7, and Cor. 2.17], we need to check that the following 
criteria are verified for the domain V of a: 

• / G V =S> / G V and o(Re/, Im/) G R, 

• fGVJ real-valued => |/| G V and a(|/|, |/|) < a(/, /), 

• < / G V =>• 1 A / G V and a(l A /, (/ - 1)+) > 0. 

It is clear that k G H 1 ^, 1) if k G H 1 ^, 1). Further, if k is real valued, then \k\ G ff^O, 1) 
and \k\' = sign*; • k' , and if < k, then 1 A k G H 1 ^, 1) with (1 A k)' = fc'l{fe<i} and 
((k-l)+y = k'l {k>1} . 

By definition, the subspace V contains exactly those functions on the network that are 
continuous in the vertices (see (|2 . 2f) ) . Take any / G V. By definition we have fj = (f)j, 
1 < j < to. It follows from the above arguments that / G (i? 1 (0, 1)) , and one can see 
that the continuity of the values attained by / in the vertices is preserved after taking the 
complex conjugate /. Hence, / e V. Moreover, a(Re/, Im/) is the sum of to integrals. 
Recall that the weights are real-valued, positive functions. Since all the integrated functions 
are real- valued, it follows that a(Re/, Img) G R. Thus, the first criterion has been checked. 

Moreover, if / is a real- valued function in V, then = \f\j, 1 < j < m, and one sees 
as above that |/| G V. In particular, ||/|'| 2 = |/'| 2 , and there holds 

m „i 

0(1/1, l/l) = E / Ho j (x)\f j (x)\ 2 dx = a(f,f). 

This shows that the second criterion applies. 
Finally, take < / G V. Then 

GMI> , 

with all the functions 1 A fj G H 1 ^, 1), hence 1 A / G (i^ 1 (0, 1))™. Again, the continuity of 
/ in the vertices imposes the same property to the function 1 A /, i.e., 1 A / G V. Further, 
there holds 

m „i 

0(lA/,(/-l)+) = ^c j (lAf j )'(x)((f j -l)+y(x)dx 

3=1 J ° 
m „i 

= X] / ^j C 3fj( X )Mf j <l}( X )fj( X ) 1 {f j >l}( X ) dx = - 
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We have checked also the third criterion, thus the claim follows. □ 

Lemma 3.2. The semigroup (T2(t)) t >o on X2 associated with a is ultracontractive. In 
particular, it satisfies the estimate 

(3.2) ||Ta(t)/|U. < Mt-i\\f\\ X2 for all t G (0, 1], / G X 2 , 

for some constant M . 

Proof. The form norm || • j| a on V is equivalent to the norm || • || y, cf. the proof of Lemma l2.ll 
Thus, by |Ou041 Thm. 6.3 and following remark] it suffices to show that there holds 

||/||x 2 < MH/lll -H/H^ for all feV, 
for some constant M. Recall the Nash inequality 

1 1 

/ 33 j II*IU 2 (0,1) ^ M l (ll fc 'llL2(04) + llfclliHO,!)) 3 ' ll fc llli(0,l) 

< ^illfc|l| 1( o,i)-!l fc !li(o,i)' 

which is valid for all k G fl" 1 (0, 1) and some constant Mi, cf. [Ma851 Thm. 1.4.8.1]. 
Take finally / G V and observe that by (|3.3p 

m rn ^ 

ll/llx 2 = ll^'Hi 2 (o,i;w^) - M ^ X ll^illffM . 1 ;^^) ' H^llii(o,i;w^) 
i=i 3=1 

2 

(m \ 3 / m 

< A/ 2 ||/||J -H/ll^, 

using the Holder inequality. Thus, the claim follows. □ 

The following now holds by |Da89| Thm. 1.4.1, Thm. 1.6.4, and Thm. 2.1.5] and |Ou04| 
Thm. 3.13]. 

Corollary 3.3. The semigroup (T2(t)) t >o extends to a family of compact, contractive, pos- 
itive one-parameter semigroups (T p (t))t>a on X p , 1 < p < 00. Such semigroups are strongly 
continuous if p G [1, 00), and analytic of angle ^ — arctan 2^/ p -i f or P ^ °°)- 

Moreover, the spectrum of A p is independent of p, where A p denotes the generator of 

(T p (t))t>Q, 1 < p < 00. 

The estimate on the analyticity angle of (T p (t)) t >o is not sharp, cf. |Mu05j for details. 
Remark 3.4. Consider the part A of A in (C[0, l]) m , whose domain is given by 

D(A) = {/ G (C 2 (0, l)) m : 4>+/'(0) = Kf'O) 

3d G C" s.t. ($+) T d = /(0) and ($") T d = /(l)} . 

Define 

C(G) := {/ G (C[0, l]) m :3rfeC" s.t. ($+) T d = /(0) and (^) T d = /(l)} , 

which can be looked at as the space of all continuous functions on the graph G. It is 
easy to see that D(A) = C{G). By Corollary 13.31 A has positive resolvent, and it follows 
by |ABHN0T1 Thm. 3.11.9] that its part in C(G) generates a positive strongly continuous 
semigroup. 
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In the next lemma we show that the generators of the semigroups in the spaces X p , 1 < 
p < oo (see Corollary [23]) have in fact the same form as in Xi, with appropriate domain. 

Lemma 3.5. For all p G [l,oo] the generator A p of the semigroup (T p (t))t>o is given by 
the operator defined in (|2.3[) with domain 

D(A P ) = {/ G riTli W 2 *(0, 1; p,dx) : $+/'(0) = and 
3d G C" s.t. ($+) T d = /(0) and ($-) T d = /(l)} . 

Proof. Let us prove the claim for p > 2. We have already remarked that X p <-> X g for all 
1 < 9 < P < oo- Moreover, it follows by the ultracontractivity of (T2(i))t>o (see Lemma 
13. 2[) that X p is invariant under (T 2 (£)) t >o for all p > 2 because ii f £ X p then / <S X 2 , and 
by 02), 

||T 2 (t)/||x, <Cr.||T 2 (t)/|| Xoo <C-Mt-i\\f\\ X2 <C-Mt-i\\f\\ Xp . 

Thus, by |EN00| Prop. II.2.3] the generator of (T p (t)) t > is the part of ^ in X p . A direct 
computations yields the claim. 

For 1 < p < 2 the claim can be proven by duality, mimicking the proof of [Mu05[ 
Lemma 4.11]. □ 

Theorem 3.6. The first order problem (|3.1j) is well-posed on X p , p G [l,oo) 7 as well as on 
C(G), i.e., for all initial data f £ X p or f € C(G) the problem p.lj) admits a unique mild 
solution that continuously depends on the initial data. 

Such a solution is bounded in the time as well as (if p > 1) in the space variables. If 
further f G X p , 1 < p < oo, and Cj G C°°[0, 1], j = 1, . . . , ra, i/ien i/ie solution u(t, ■) is 0/ 
class (C°°[0, l]) m /or t > 0, and in particular the problem is solved pointwise for t > 0. 

Proof. The well-posedncss and boundedness results follow from the fact that the operators 
A p generate ultracontractive analytic semigroups. If Cj G C°°[0, 1], j = 1, . . . , m, then we 
can show as in the proof of Theorem [23] that D(A™) C (C°°[0, l]) m for all p G [1, 00]. Since 
the semigroup (T p (t))t>o is analytic, 1 < p < 00, it maps X p into D(A p °), and the claim 
follows. □ 



4 A characteristic equation 

Having proved that the Cauchy problem (|3.ip is well-posed, we want to study the qualitative 
behavior of its solutions. To this end we investigate the spectrum of the generator (A, D(A)). 
Since by Corollary 13.31 the spectra of all A p on X p , 1 < p < 00, coincide, it suffices to study 
the operator A = Ai on Xi. Hence, we are interested in the spectrum of the operator 




with domain 

(4 2) D(A) ■■= {/ G 117= 1 (H 2 (0, ly^jdx) : $.+ /' (0) = and 

3d G C" s.t. ($+) T d = /(0) and ($-) T d = /(l)} . 

Recalling properties of (A,D(A)) already yields some information on its spectrum. 

Lemma 4.1. The spectrum of (A, D(A)) lies on the negative real line and consists of eigen- 
values only. Moreover, s(A) = G o~(A). 
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Proof. First note that 1 G and Al = 0, thus A is not invertible and € cr(A). Since 

A generates a contractive semigroup (cf. Proposition 12 .4[) . s(A) = 0. It follows from |EN0Q[ 
11.4.30.(4)] that, since D(A) is contained in (i? 2 (0, 1); /ijdx) , the resolvent of A is compact. 
Therefore the operator A only has point spectrum. Recall that by Proposition 12.41 the 
operator A is self-adjoint, hence all its eigenvalues are real. □ 

From now on we will assume that all the weights Cj, i = 1, . . . , m, are constant. Our aim 
is to find a 'characteristic equation' for the spectrum of A. In particular, we will be able to 
connect the eigenvalues of the operator A to the eigenvalues of the Laplacian or admittance 
matrix of the corresponding graph. This is the n x n matrix 

(4.3) C := D — A, 

where A is the standard 0—1 adjacency matrix of the graph and D the diagonal matrix 
of vertex degrees. It is well known that its spectrum reveals many properties of the graph, 
hence it is used in many applications (see e.g. |Ch97l lMe94[ IMo97] ) . 

We further define the generalized weighted adjacency matrix of the graph G in the case 
< A 7^ CjPir 2 , j = 1, . . . , m, I e Z, as 



(Ac(A))i 




ifr( Vi )nr(v fc ) = { 
if jer( Vi )nr( Vfc ) 




By Dc(X) we denote the n x n diagonal matrix (again for < A ^ Cjl 2 n 2 , j = 1, . . . , m, I € 
Z) defined as 

D c {\) : diag 

Finally, we define the generalized weighted Laplacian matrix as 

C C (X) :=D C (A)-A C (A). 

We will now express the above matrices using the weighted incidence matrices. For this 
purpose we define diagonal matrices 

( x 

Sin x := diag sin — — , . . . , sin 



Cosx := diag ( cos — ^= cos 



Cotx := Sin • Cosx, and 
C := diag(l/VcI,...,l/ A /c^). 

Lemma 4.2. For < A ^ c r l 2 Tr 2 , r = 1, . . . , m, I € Z we have 

Ac(X) = ■C-Sin- 1 VA-($-) T + $--C-Sin- 1 VA-($ + ) T and 
D C {X) = $+ •C-Cot\/A-($ + ) T + $.--C , -CotVA-($-) T . 



We are now able to describe the spectrum of our operator A in terms of spectral 
values of £c(A). Similar results have already been obtained in much the same way as 
in [Be85] . fNlgg] . jNi87] . |Ni87bj . |Be88bj . and [Ca97] for the cases fi 3 = 1 and/or c 3 = 1. 

Theorem 4.3. For the spectrum of the operator (A, D(A)), defined in (|4.1[) - (|4.2p , we obtain 

a(A) = {0}Ua c U(T£, 
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where 



a c C {-CjfcV : k E 1\{0},i = l,...,m} and 
(T £ = {-A G R_ : A 7^ c.fcV, det£ c (A) = 0} . 



Furthermore, 



1. A = G o~(A) is always an eigenvalue of (geometric and algebraic) multiplicity 1 with 
an eigenvector f(x) = 1, i/ie constant 1 function. 

2. -A G (T£ is an eigenvalue of A with corresponding eigenvector 

f{x) = Cos^/Xx- ($ + ) T d + Sin- 1 VA-Sin\/Aa;- (($~) T - CosVA- ($ + ) T ) d 



where d G ker£c(A), so £/ie multiplicity m(— A) o/ </izs eigenvalue is equal to 
dimker£c(A); 

5. — Cik 2 ir 2 G <7c is an eigenvalue of A if and only if there exist b G C m and d G C" smc/i 
i/iai whenever j G r(v r ) n T(v s ), j G {1, . . . , m} 7 we /lave 



(4.4) 



d r = (-l)^? fe d s , if^keZ, 

bj — sraT 1 . / ^-fcTr • d r — cot -1 , ^kir ■ d s , otherwise. 



These vectors further satisfy the equation 

(4.5) • C ■ Sin y^wr • ($ + ) T d = ■ C ■ Cos ^kir - ■ C) ■ 6. 

// t/ie eigenvector exists, then it has the form 

f(x) = Cos yfcik-Kx ■ ($ + ) T d + Sin ^/clknx ■ b. 

Proof. By Lemma r4.il we need to solve the equation 

Af = -\f for / G D(A) and A > 0. 

We will distinguish three cases. 
Case 1: Assume that A ^ Cik 2/ n 2 for all k G Z, i = 1, . . . , m. 
In this case the eigenfunctions of A are of the form 

f(x) = Cos V~\x ■ a + Sin \T\x ■ b for some a, b G C"\ 

From the continuity assumption in the domain of A (see (|4.2[) ) 

3d G C n s.t. ($+) T d = /(0) and ($~) T d = /(l) 

we obtain 

f(x) = CosVAx- ($+) T d + Sin- 1 \/A- SinVAx- (($~) T - CosVA- ($+) n 

for some d G C™. The other condition $+/'(0) = in the domain £>(A) (i.e. the 

Kirchhoff law) yields that / G ker (A — A) if and only if the vector d G C™ satisfies 



$+ • C- Sin^VA- -CosVX- ($+) Jd = 

= C - (CotVA- (($-) T - CosVA - ($ + ) T ) -Sin\/A- ($+) t ) d. 
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Observe now that, by Lemma [4.2[ the following two terms are the previously defined diagonal 
matrix 

• G ■ Cot VA ■ ($+) T + • C ■ Cot\/A ■ ($") T = D C (X), 
while rearranging the remaining terms yields the weighted adjacency matrix 

$+-C-Sm- 1 V\-($-) T + &--C- (Sin^x/A- Cos 2 \/A + Sin \/a) • ($+) T = 
= *+ • C • Sin ^VX ■ ($~) T + • C* • Sin • ($+) T = A C (A). 

Summing up, the condition for d G C" becomes 

(A c (A) - £>c(A)) d = 0, that is d G kcr£c(A). 

Case 2: A = 0. 

The eigenfunctions of A, corresponding to A = 0, are of the form 

f(x) = x ■ a + b for some a, b G C m . 

We repeat the above procedure and the conditions in the domain (|2.4[) yield 

f(x) = ($+) T d - x • $ T d for d G ker$„,$ T with = - 

Since our graph is connected, the multiplicity of in a ($„,$ T ) is 1 (cf. [Ch97l Lemma 
1.7. (iv)] or |Mo971 Prop. 2.3])). It is easy to see that the corresponding eigenvector equals 
d = 1 := (1, . . . , 1) T . Now compute 

f(x) = ($+) T l - x • $ T 1 = 1 for all x. 
Case 3: A = Cik 2 ir 2 for some k G Z \ {0} and some i G {1, . . . , m}. 

We proceed as before while some care need to be taken with zero entries that arise. Before 
applying the inverse of Sin y/\, the continuity condition in (|4.2|) implies 

f(x) = Cos y^Cikirx ■ ($ + ) T <i + Sin yfcikitx ■ b, 

where b satisfies the equation 

(4.6) Sin VcI/ot • b = ((<J>~) T - Cos^ctf • ($ + ) T ) d. 

Since the i-th entry on the left-hand side equals 0, the vector d should satisfy the condition 

d. p = (-i) k d q for i g r(vp) n r(v 9 ). 

Moreover, if for any other j G {1, . . . , m\ we have ,/—k G Z, then also 

y c j 

d r = {-i)y^ k d s for j g r(v r ) n r(v s ). 

For each of these cases we have no conditions on bj . If on the other hand 4- Z, equation 

(gH) yields 

bj = sin -1 , —kn ■ d r — cot -1 . —kn ■ d S7 j G r(v r ) l~l T(v s ). 
V c i V c 3 

Furthermore, the condition $+/'(0) = <F~/'(1) in the domain D(A) (i.e. the Kirchhoff law) 
implies that above vectors d and b have to satisfy the equation (|4.5I) . □ 



12 



Let us emphasize that the condition 3 in the above theorem in not always satisfied, 
therefore the spectrum of our operator A strongly relies on the underlying graph and on the 
weights Cj. 

From now on we assume cy = 1, j = 1, . . . , m. In this case we are able to connect the 
spectrum of the operator A to the spectrum of yet another matrix known in graph theory 
The transition matrix is defined as 

P := D- X A 

and is studied in connections with random walks on graphs (see e.g. |Mo97[ Sec. 5.2] or |Ch97| 
Sec. 1.5]). The matrix P is always positive, symmetric, row stochastic matrix with eigenval- 
ues 

a (P) = {ai, . . . , a„} where — 1 < a n < • • • < ct2 < ol\ = 1. 

By |Mo97| Claim 5.3], 1 is a simple eigenvalue whenever G is connected (what we assumed 
at the beginning) and — 1 is an eigenvalue if and only if G is bipartite (see also Lemma 1 in 
Section 5 of |Be85 ). It turns out that an important subset of the spectrum depends on the 
fact whether the graph G is bipartite or not. This property means that the set of vertices 
can be divided into two disjoint subsets V\ and V2 such that any edge of G has one endpoint 
in one and the other endpoint in the other subset. Note that G is bipartite if and only if it 
does not have any odd cycle. 

The following characteristic equation has already been proved by von Below jBe85| . We 
state and sketch the proof in our context for the convenience of the reader. 

Theorem 4.4. Let (A, D(A)) be the operator defined in ()4.1|) - (|4.2|) . with Cj = 1, j = 
1, . . . , m. Then for the spectrum of A we have 

a(A) = {0}Ua p Ua k , 

where 

a p = j- (22tt± arccosa) 2 : a G <r(P) \ {-1,1} and I G zj , 

and 

cr fe = {-fcV :fceZ\{0}} 
For the multiplicities of the eigenvalues we have: 

1. m(0) = I; 

2. m(— A) = dimker ^cos-y/A ■ 7 — PJ for —A 6 a p : 

3. m(—k 2 ir 2 ) = m — n + 2, if G is bipartite; 

4- m(— 4l 2 ir 2 ) = m — n + 2 and m(—(2l + l) 2 7r 2 ) = m — n, if G is not bipartite. 

Proof. We use Theorem 14.31 for C = I. For the spectral point A = the statement follows 
directly from Theorem 14.31 Assume first that A ^ k 2 ir 2 for any k e Z. Then 

£r(A) =sin~ 1 \/A(cosv / A-.D-A) , 

and the characteristic equation becomes 

-A e cr (A) <^=S> det (cos VA • D - A\ =0 <*=^> dot (cos V\ ■ I - p) = 

for the transition matrix D~ l A = P. The last equivalence says that 

—A G a (A) cos VA G a (P) A = (21 7r ± arc cos a) 

for some a G a (P), — 1 < a < 1, and / G Z. Since dimker(cosv / A--C' — A) = dimker(cos V\ ■ 
I — P), statement 2 also follows by Theorem 14.31 2. 
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Now, let A = fc 2 7r 2 for some k £ Z \ {0}. Observe that condition (|4.4p in Theorem 14.31 3 
becomes 

d r = (-l) k d s whenever T(v r ) n r(v s ) ^ 0. 

If k is even this condition is always fulfilled for d = c ■ 1 for any cgl. Because the network 
is assumed to be connected, there is no other solution. For odd k we can always choose 
d = 0. However, we can find a nonzero d only in the case G does not have any odd cycles, 
that is when G is bipartite - hence, when his set of vertices can be divided into two disjoint 
subsets Vi and Vi such that any edge of G has one endpoint in one and the other endpoint 
in the other subset. If this holds, the coordinates of d can be chosen in such a way that 
at the places of vertex indices belonging to V\ we set c and at the places of vertex indices 
belonging to V 2 we set — c, c£l. By connectivity this are again all possible solutions. 
Since sin kir = and cos kn = (— l) fc , the other condition (|4.5p in Theorem l4.3l 3 becomes 

($i-(-l) fc $-)6 = 0. 

Using the proof of |Be851 Theorem 5 (17)] we obtain 3. and 4. □ 

5 Stability results for the diffusion problem 

In the last section we are interested in the asymptotic behavior of solutions to the problem 
(|3.1[) . By Corollary 13. 3[ the corresponding semigroup T p (t) t >o on X p = (L p (0, 1); fijdx) m , 
1 < p < oo, has many nice properties: it is contractive, compact, positive. These properties 
already yield norm convergence of the solutions to an equilibrium (cf. |EN00( Cor. V.2.15]). 

From the connectedness of our graph, used in the proof of Thm. 4.3.1, we obtain another 
useful property of the semigroup. 

Proposition 5.1. The semigroup (T p (t)) t >Q on X p , p £ [l,oo) ; is irreducible. 

Proof. It is enough to prove the statement for p = 2, because the irreducibility is inherited 
for the extrapolation semigroups in X p , 1 < p < oo, using Corollarv l3.3l and [Ar041 Theorem 
7.2.2]. Since X 2 is reflexive and (T 2 (t)) t>0 is bounded, by jENOOl Example V.4.7] we have 
that the semigroup is mean ergodic (see |EN0Q[ Definition V.4.3]). By Theorem 14.31 1. we 
obtain that the corresponding mean ergodic projection P defined by 

1 f 

Px := lim — / T(s)xds 

r^oo r J 

is the projection 1 (g) 1 onto the subspace <1>. Let {0} ^ Jc X 2 be a closed invariant ideal 
for (T(t)), that is T(t)J C J, t > 0. Then also PJ C J holds. By definition, PJ C<1> and 
so J contains the closed ideal generated by a constant function - hence the whole space X 2 . 
From this follows that the semigroup is irreducible. □ 

Knowing that our semigroup is irreducible we can now show its norm convergence towards 
a projection of rank one. 

Corollary 5.2. For the semigroup (7 1 p (<)) t > on X p , p £ [l,oo), the following hold. 

1. The limit Pf := limt_>. 00 T p (t)f exists for every f £ X p . 

2. P is a strictly positive projection onto <1>= KcrA, the one- dimensional subspace of 
X p spanned by the constant function 1 . 

3. For every e > there exists M > such that 

(5.1) ||T p (£) - P\\ < Me [e+X2)t for all t > 0, 

where X 2 is the largest nonzero eigenvalue of the generator A. 
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Proof. The first assertions follows from |Na861 Cor.C-IV.2.10] and the second one from 
|Na86[ C-III.3.5.(d)] and Theorem 14.31 Since P is the residue corresponding to the spectral 
value Ai = and is a first order pole of the resolvent, estimate (3.2) in [ENOCH Cor. V. 3. 2] 
yields the third assertion. □ 

Note that the property of converging to an equilibrium does not depend on the structure 
of the network. However, the speed of the convergence towards a projection is determined 
by the second largest eigenvalue A2 of A and thus relies on the network. 

Combining graph theory and results about the spectrum of A, obtained in Section 01 
we can now draw some further estimates containing graph parameters for the speed of 
convergence in (|5.1I) . Let us demonstrate this for the case when all Cj = 1 and the graph is 
regular, that is, every vertex has the same degree. This means that 

|r(vi) I = 7 for alH = 1, . . . , n, and D = 7 • I. 

Two generic examples of regular graphs are the n-cycle C n and the complete graph K n on 
n vertices (in the latter, every two vertices are connected by an edge). 
In this case the characteristic equation becomes 

-Ae<r (A) det (cos \/A • 7 • / - A ) =0 cos VA • 7 e a (A) , 



see the proof of Theorem 14.41 or [Bc85, Sec. 6]. Using the Laplacian of the graph, defined in 
(|4.3[) . we obtain for its spectrum 

v&a(£.) det (v ■ I - (7 • I - A)) = <^=> -v + 7£ cr(A). 

Hence, investigating the spectrum of the generator A, we are looking for A's such that 

\ \ 2 

v 



(5.2) A = - ^2/7r±arccos \\ - -J J , ZeZ 

where v G c(£). The spectrum of C is sometimes also called the spectrum of the graph 
and is well investigated in graph theory In the following we always refer to the survey 
paper |Mo91j . 

Example 5.3. In the case when G = C n , the eigenvalues of C are precisely the numbers 

„ „ /2(fc-lW\ , , 
Vk = 2 — 2 cos , k = 1, . . . , n, 



while for G = K n we have 

v\ = and Vk = n for 2 < k < n. 

As explained above we are interested in the second-smallest eigenvalue V2 € o~ (£), which 
is also called the algebraic connectivity of the graph, see |Mo91j . It is related to the classical 
connectivity parameters of the graph (see below). If we look at (|5.2p we can easily conclude 
that 

(5.3) A2 = — ^arccos ^1 - 

since the function arc cos is strictly monotone decreasing and assumes its values between 
and 7r. 
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Example 5.4. By increasing the number of vertices n, the convergence gets slower on the 
cycle C n and on the complete graph K n . In fact, for C n we have 7 = 2, hence by Example 
15.31 we obtain that 

f . 2tt,\\ 2 An 2 

A? = — arc cos 1 — ( 1 — cos — 



nil n 2 



For K n we have y = n — 1, hence by Examvle \5.3l 

At = — ( arc cos ( 1 ) ) = — ( arc cos ( — 



n — 1 J J \ \ n — 1 

If we have an estimate from below for u%, using (|5.1|) and formula (|5.3p . we obtain an 
upper estimate for the convergence speed of the semigroup to the one-dimensional projection. 
In |Mo91j we find many estimates for v 2 that use several graph parameters. As an example 
we mention the next. 

Definition 5.5. The edge connectivity parameter n = 77(G) of the graph G is defined as 
the minimum number of edges whose deletion from G disconnects it. 

Example 5.6. By \Mo91[ Theorem 6.2(b)], 

7T 

v 1 > 2*7(1 — cos — ) 
n 

holds, where n is the edge connectivity parameter of the graph. From this we obtain for (|5.ip 
that for every e > there exists an M > such that 

\\T P (t) - P\\ < Me (£ - (arccos(1 -^ (1 - cos " )))2)t for all t > 0. 
Another estimate can be obtained using the diameter of G. 

Definition 5.7. The diameter of G denoted by diam(G) is the largest number of vertices 
which must be traversed in order to travel from one vertex to another when paths which 
backtrack, detour, or loop are excluded from consideration. 

Example 5.8. By \Mo91[ (6.10)], we have 

4 

u 2 > 



n ■ diam(G) 

Hence, again we can conclude that for every e > there exists an M > such that 
\\T p (t) - P\\ < Me (£ ~ (arccos(1 ~ ^" ^<°> ))2)t for all t > 0. 
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